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Direct constraints on charged excitations of dark matter 
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If the neutral component of weak-scale dark matter is accompanied by a charged excitation 
separated by a mass gap of less than ~20 MeV, WIMPs can form stable bound states with nuclei. 
We show that the recent progress in experiments searching for neutrinoless double-beta decay sets 
the first direct constraint on the exoergic reaction of WIMP-nucleus bound state formation. We 
calculate the rate for such process in representative models and show that the double-beta decay 
experiments provide unique sensitivity to a large fraction of parameter space of the WIMP doublet 
model, complementary to constraints imposed by cosmology and direct collider searches. 



Introduction Weakly interacting massive particles 
(WIMPs) are well-motivated candidates for dark matter 
(DM), that offer a variety of potential non-gravitational 
signatures, including the possibility for a laboratory de- 
tection. Model realizations of WIMP DM often entail 
multiple states. Frequently encountered are scenarios in 
which DM (X°) is part of a multiplet with an electri- 
cally charged excited state (X^) that enables to natu- 
rally regulate its abundance through co-annihilation. As 
has been noted in [l| , if the mass gap between the neutral 
and charged state is sufficiently small, MeV-scale energy 
depositions become possible via the process of WIMP- 
nucleus bound state formation. Such models are of spe- 
cial interest because they also lead to an enhancement 
of WIMP annihilation to mono-energetic 7-rays [l| , and 
well-defined collider signatures through the production of 
the charged WIMP states. 

Denote by (NX~) the bound state of DM with a target 
nucleus N. Depending on the relation between spins of 
X and X~ two generic scenarios can be envisaged, 



Case A: N z 
Case B: N z 



X° 
X° 



(N Z X-) 
(Nz+iX- 



(1) 
(2) 



where Z denotes the charge of N. If (NX~) is not in its 
ground state, it will de-excite by emitting 7-rays. Thus, 
the observables for these processes are the positron and 
7-rays (and a residual anomalously heavy nucleus.) The 
recombination is guaranteed to happen once the (NX~) 
Coulomb binding energy Eb allows to bridge the mass 
gap Am = m x - —m x o in the DM multiplet. Mass split- 
tings of ^20 MeV or less can be probed through (TTJ) and 
P|) . Both scenarios are readily realized in supersymmet- 
ric (SUSY) scenarios. 

The primary purpose of this letter is to show that the 
recent advances in neutrinoless double-beta (0^/3/3) decay 
experiments allow for a direct probe the neutral-charged 
WIMP doublet models. The 0(1-10 MeV) scale energy 
deposition in ((T|) and © give Qv(3(3 searches a clear ad- 
vantage over direct DM detection experiments. The lat- 
ter are optimized to detect keV-scale energy depositions 
from the elastic scattering of WIMPs and the MeV-scale 
energy range is usually not reported. Moreover, albeit 



similar detection principles and background rates, CV/3/3 
experiments are more readily scaleable in mass, and the 
next decade is bound to produce progress not only in 
limiting/detecting the Majorana mass for neutrinos (the 
primary goal of Ov/3/3 programs), but also in constrain- 
ing the charged excitations of WIMPs. In what follows 
we consider representative models for both cases and use 
them to set constraints from new data released this year 
by EXO-200 Q and Kamland-Zen Q Ov/3/3 collabora- 
tions. 

Representative Models Case A demands that the 
spins in the DM multiplet differ. With an eye on SUSY 
we take the DM particle to be a Majorana fermion 
X° = x (like the neutralino) and the charged excited 
state to be a complex scalar X~ = f~ = f (like the stau.) 
The effective Lagrangian which governs the bound state 
formation is of Yukawa type 



C A = x(9clPl + 9eR^R)ef f + h.c. 



(3) 



with in general complex couplings g e L,R for the chirality 
projections Pl,r- 

In Case B the DM states share the same spin. For sim- 
plicity, we take both as scalars in the electroweak doublet, 
a real scalar X° = v° (like a sneutrino) and a complex 
scalar X~ = f~ as before. The relevant effective La- 
grangian can be written as 

C B = ?fw->*(d^i> -^d li P )+h.c. + C coat .(A) 

where £ C ont. contains contact interactions; in SUSY 
ffcff = 92 cos Of cos Qya , where (72 is the gauge coupling 





(a) (b) 
FIG. 1: Illustrations of case A (a) and case B (b). 
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TABLE I: Relevant nuclei, binding energies, and expo- 
sures MT for different experiments considered in this work. 



Exps. 


EXO-200/ 
Kami. -Zen/ 
XelOO 


DAMA 
Nal(Tl) 


SNO 
NaCl 


Bore- 
xino 


Nucleus 


Xe 


I 


CI 


C 


Ei 0) (MeV) 


18.4 


18.2 


6.3 


2.7 


MT (kg yr) 


40/30/0.9 


7.5 


1274 


1.3 x 10 5 



of the weak interaction, Of is the mixing angle between 
the left and right handed stau, and 9„o is the mixing an- 
gle between the sneutrino and any-sterile sneutrinos in 
the theory. 

Bound state properties and formation In Tab. HI we 
list the ground state energies for the X- bound 

states with relevant nuclei, using the homogeneous nu- 
clear charge distribution within sphere of radius Rq = 
x/5/3R vms , and i? rms values are taken from Ref. Ki- 
netic energies 0(100 keV) of the incoming DM and that 
of the recoiling bound state can be neglected, and the ex- 
pected total visible energy Etot injected in the detector 
is given by 




Case A, 
mz — mz+i Case B. 



(5) 



(n,l) 



In Case A the positron is emitted with an energy E, 

_ _ m ^ w here m e is the electron mass and n, I 
denote the usual initial principal and orbital quantum 
numbers of the capture. After the positron is stopped, its 
annihilation with an electron yields an additional energy 
of 2m e whereas the excited (Nf)* relaxes by emission of 
gamma-rays. In case B no positron is emitted but the dif- 
ference in nuclear mass upon the nuclear transmutation 
Nz — > Nz+i becomes accessible. 

For case A, the product of recombination cross sec- 
tion a a with the incoming DM velocity v reads, 



a A v 



(1.9, 



eL 



\geR\ 2 )/(8Trm x )xj:.B n ^ (6) 



where B n j is the contribution from capture into state 
(n, I). In the limit Am 3> m e and p x 3> p e + , 



B n ,i ~ (e^ i1) - Am -m e j y v ^ b 

x [ d 3 r 1 d\ 2 r n i(ri)<f>nAr2)e^< r ^ ) . (7) 



Here, is the wave function of the relative motion of 
(Nf) with reduced mass fi. For n < 50 we calculate B n ^ 
explicitly by numerical solution of the Schrodinger equa- 
tion. Whenever Am is small, capture into a multitude of 
states (also with n > 50) is possible. Importantly, in this 



case we expect a a to flow towards a semi-classical (SC) 
limit: when \ approaches the nucleus, a critical radius 
T6 in the electromagnetic potential V(r) will be reached 
when the transition x ~ ► t~ + e + becomes energetically 
possible, V(rb) + Am + m e < 0. The integration over the 
fly-by time when r < r& then gives the rate for this tran- 
sition to happen. In this SC limit we find for ^ n , B n j, 



r\<r b 



d 3 fx/(V(\r\) + Am) 2 - m 2 e (V(\r\) + Am ■ 



m e ). 
(8) 

Whenever, say, n > 10 is accessible in the capture, we 
find a perfect agreement between SC and the explicit 
quantum mechanical calculation. For lighter nuclei with 
shallow binding the latter is the preferred choice as typ- 
ically only few states contribute to the signal above the 
detector threshold. As expected, being an inelastic pro- 
cess, ctav is largely independent on v, and to very good 
approximation (ctav) ~ o~av. 

Once Am and DM mass m x are chosen, the induced 
signal in the neutrino searches can be translated into 
a constraint on the combination of Yukawa couplings 
(ISeil 2 + \gen\ 2 )- We choose to trade the latter against a 
constraint on the f lifetime instead, or, alternatively on 
its decay width F? = Tf . In the same approximation as 
above, 



47T77K 



-(Am - 



eL 



\9eR\ 2 ). (9) 



Let us now consider case B. During the capture an in- 
terconversion from a neutron (n) to a proton (p) takes 
place inside the nucleus. Without detailed knowledge of 
the nuclear wave functions, an exact QM calculation is 
not possible, although it is expected that the cross section 
can be dominated by a series of narrow resonances espe- 
cially when Am is large [3, [BJ . When many bound states 
are available for capture, one expects the SC method to 
give a statisfactory description, which we employ below. 

We use the Fermi gas model for the density and mo- 
mentum distributions of n and p inside the nucleus Nz- 
The calculation of o~b proceeds in two steps. First, we 
compute the fundamental process of n — > p in the pres- 
ence of v° ' . Because of Pauli blocking, part of the binding 
energy is invested in elevating p above the Fermi surface 
of Nz+i- In a second step we again obtain the total cross 
section by consideration of the fly-by time of u° in which 
the transition is possible. The result reads, 



ffcff m p 

8M^ 



d 3 rp n (r) 



dpi 



\J (-V(r) - Am + p° ) 2 - m 2 p \Jp° n 2 - m 2 n 
xe(-V(r)-Am +P :-m p - P ^[ Nz+iy 



(10) 
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where p n is the number density of neutron inside the 
target nucleus and p p f(Nz+i) is the Fermi momentum 
of the Nz+i nucleus. The latter is determined from m p + 

- E£ 75 = m n + ! - a ^ 1 — -. tor the capture on xenon, 

N z = Xe and N z +\ = Cs. 

Neutrino experiments We now come to the central 
part of our study — demonstrating the potential of neu- 
trino experiments and in particular of 0^/3/3 searches in 
probing DM. As we have argued, the total energy in ([5]) 
is monochromatic with all its energy injected almost in- 
stantaneously. The signal shape is hence determined by 
the energy resolution o~ e of the experiment. Given an ex- 
posure MT the total number of events in an energy bin 
AE = E max - E min then reads, 



MTN T pD M (o-v) 



Erf 



2mxo 

^max -^tot 

V2a e 



-Erf 



-^min — E t , 

V2a e 



(11) 



where pum ~ 0.3 GeV/cm 3 is the local DM energy den- 
sity, Nt is the number of target nuclei per kg of active 
detector mass. 

We start by considering the EXO-200 experiment Q 
which has 110 kg fiducial mass enriched in 136 Xe, whereas 
for our purpose the entire mass of Xe in the fiducial vol- 
ume becomes active, and the exposure as a DM target 
is listed in Tab. H The TPC-type setup of EXO-200 al- 
lows to distinguish between events occurring at a single 
sites (SS) or at multiple sites (MS). The capture typically 
qualifies as an MS event, given the macroscopic range of 
produced 7-rays (~ 6 cm at 1 MeV in LX.) We use the 
last bin containing 24 events above 3.5 MeV to set a con- 
straint which thereby probes the region Am < 15 MeV. 
Two complications arise: 1) part of E to t can be taken 
out of the fiducial volume by 7-rays and 2) events at SS 
with depositions greater than 10 MeV are automatically 
discarded. We have studied both effects in a dedicated 
Monte Carlo (MC) analysis and find that this can weaken 
the limits by up to a factor of two and which is taken into 
account in presented constraints; details will be presented 
in The results for cases A and B are shown in Fig .H\ 

The next experiment we consider is Kamland-Zen [3j 
which consists of a LX balloon submersed into the Kam- 
land neutrino detector. To obtain a conservative con- 
straint, we use the reported spectrum from which we only 
subtract the background from the 2/3 decay of 136 Xe. The 
energy resolution is excellent so that we only need to con- 
sider the bin in which E tot lies including the neighboring 
ones (with appropriate statistical penalty) in deriving the 
constraint. Unfortunately, Kamland-Zen docs not report 
a cumulative last bin like EXO-200 does. Therefore, we 
can only probe energy depositions below the last reported 
bin at 3.8 MeV. This, in turn, restricts the considered 
mass splittings to Am > 14 MeV as shown in Fig. [3] 



Once Am becomes large enough to approach the thresh- 
old, the limits become correspondingly weaker because of 
limited phase space and higher count rates. 

Turning to small values Am < 4 MeV we either probe 
energy depositions beyond 10 MeV with heavy elements 
like with Xe or, alternatively, we can probe smaller (and 
hence reported) energy depositions using lighter nuclei. 
For the latter, the neutrino experiments of choice are cur- 
rently the salt phase of the SNO water Cherenkov detec- 
tor with a 2 ton loading of NaCl and the Borexino detec- 
tor with its carbon-based scintillator. We set constraints 
using the respective event spectra reported in [§] and [13] 
where in the former case we subtract the solar neutrino 
signal; see Q . Finally, we note that Am < 4 MeV is also 
challenged by severe limits [ll[ of anomalously heavy C, 
N, O nuclei among others. From Fig. [2] it can be seen 
that the limit from Borexino is in the same ballpark. 

Direct detection experiments One of the current most 
sensitive direct detection experiments is XENON100. 
Since it is a LX detector, we can directly compare it 
to EXO-200 and Kamland-Zen. Using the electromag- 
netic background spectrum up to 3 MeV published from 
a commissioning run 12j, we find that constraint is at 



least one order of magnitude weaker. This can be traced 
back to higher levels of radioactive background present 
in XENON100. 

A competitive experiment in sensitivity can be 
DAM A/LIBRA with its 250 kg of Nal scintillating crys- 
tals. A high-energy "single-hit" spectrum from 9 out of 
25 crystals of a 24 day run was reported in [13j . No events 
above 10 MeV were observed. A conservative constraint 
in case A is obtained by considering only recombination 
into the IS ground state of I. No 7-rays are produced 
which may well escape one crystal challenging a "single- 
hit" observation. For this reason we refrain from studying 
case B which certainly requires MC modeling. For case A 
we find that the constraint is superseded by the one of 
EXO-200. We stress that if DAM A were to present us 
with a high-energy spectrum of a larger exposure, a very 
strong result could be obtained. 

Constraints from Colliders For heavy charged sta- 
ble particles, the CMS group placed a very stringent 
bound based on the data with an integrated luminos- 
ity of 5 fb _1 and a centre-of-mass energy of 7 TeV 14 1. 



Events were collected based on two triggers. For the first 
one, the heavy stable candidate events were defined as 
tracks reconstructed in the inner tracker detector with 
large dE/dx and high pt, where dE/dx is the rate of en- 
ergy loss due to the ionization. For the second one, the 
tracks had to be additionally identified as muons by the 
outer muon detectors, thus putting a lower bound on the 
time-of-flight. 

For case A we obtain a limit on F? from the constraint 
on the LHC f production cross section <7 pro duction, 



^productions v -r~i v T < ^constraint 



(12) 
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FIG. 2: Left: Constraints for "Case A" from underground rare event searches as labeled. Darker (lighter) shadings are for 
m x = 1 TeV (100 GeV). Also shown are lower limits from CMS searches on long-lived charged states. Right: Corresponding 
constraints for "Case B". The dotted lines are contours of constant f lifetime; lighter shadings are for m x = 200 GeV. 



where R is the transverse distance, vt the transverse ve- 
locity and 7 the boost factor; ^constraint we take from [l4j |. 
The lower limits on IV for ra x = 100 GeV (and above) 
are shown in the left Fig. [5J The combined constraints 
rule out the entire region Am < 5 MeV as long as 
m x < 210 GeV; larger values of m x are currently not 
constraint by CMS. In the case that the excited state is 
a spinor, the constraint will become stronger due to that 
the dominant production channel is s-wave, instead of p- 
wave in the scalar case. ATLAS also conducted searches 
for the quasi-stable charged particles. Though similar 
sensitivity can be expected, the reported results are spe- 
cialized to a particular SUSY-breaking framework, and 
we do not use it here. In case B, f escapes the detector 
because T? is suppressed by Am 5 /m^. The CMS con- 
straint is hence trivial: rap > 210 GeV. Finally, for life- 
times Tf > 10 3 s a primordial abundance of f ~ (which is 
expected to be similar to that of v° since Am /rap <1%) 
will have a large effect on the primordial light elements; 
see [H| and references therein. The approximate con- 
straint is labeled "CBBN" . 

Outlook We have shown that rare underground event 
searches can be used to obtain strong limits on charged 
excited states of DM with mass splittings Am < 20 MeV. 
Significant improvements on these bounds can be ex- 
pected from nearly all upcoming 0v>(3/3 searches such 
as EXO, GERDA, CUORE, MAJORANA, NEXT, Su- 
perNEMO among others (see [l6| and references therein), 
as well as from future data sets of Kamland-Zen and 
EXO-200. In addition, the Gd- loadings in Super- 
Kamiokande and SNO+ will provide sensitive probes in 
the intermediate Am regime. In some cases great sensi- 
tivity is already experimentally established but unfortu- 



nately not published. For example, Kamland-Zen has not 
reported its count rate above 3.8 MeV. Presumably, the 
associated constraint can be the dominant one through- 
out the entire considered region in Am. Likewise, DAMA 
has shown a fraction of its high energy spectrum with an 
exposure of 24 days — whereas more than 100 times of this 
is available. We urge these collaborations to present this 
data. A lot of progress can be expected in the coming 
years. 
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